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IS STRONG SASI ACTIVITY THE KEY TO SUCCESSFUL NEUTRINO-DRIVEN SUPERNOVA EXPLOSIONS? 
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ABSTRACT 

Following a simulation approach of recent publications we explore the viability of the neutrino-heating explo- 
sion mechanism in dependence on the spatial dimension. Our results disagree with previous findings. While we 
also observe that two-dimensional (2D) models explode for lower driving neutrino luminosity than spherically 
symmetric (ID) models, we do not find that explosions in 3D occur easier and earlier than in 2D. Moreover, we 
find that the average entropy of matter in the gain layer hardly depends on the dimension and thus is no good 
diagnostic quantity for the readiness to explode. Instead, mass, integrated entropy, total neutrino-heating rate, 
and nonradial kinetic energy in the gain layer are higher when models are closer to explosion. Coherent, large- 
scale mass motions as typically associated with the standing accretion-shock instability (SASI) are observed to 
be supportive for explosions because they drive strong shock expansion and thus enlarge the gain layer. While 
2D models with better angular resolution explode clearly more easily, the opposite trend is seen in 3D. We 
interpret this as a consequence of the turbulent energy cascade, which transports energy from small to large 
spatial scales in 2D, thus fostering SASI activity. In contrast, the energy flow in 3D is in the opposite direction, 
feeding fragmentation and vortex motions on smaller scales and thus making the 3D evolution with finer grid 
resolution more similar to ID. More favorable conditions for explosions in 3D may therefore be tightly linked 
to efficient growth of low-order SASI modes including nonaxisymmetric ones. 
Subject headings: supernovae: general — hydrodynamics — stars: interiors — neutrino 



1. INTRODUCTION 

Recent simulations in two dimensions (2D) with sophisti- 
cated neutrino transport have demonstrated that the neutrino- 
driven mechanism, supported by hydrodynamic instabilities 
in the postshock layer, can yield supernova explosions at 
least for some prog enitor stars (11.2 and 15M G ones in 
iMarek & Jankdl2009llMuller et al.l)2012l) . The explosions oc- 
cur relatively late after bounce and tend to have fairly low 
energy, being "marginal" or only slightly abov e the "criti- 
cal threshold" in this sense. ISuwa et aTT d2010t) obtained a 
similar explosion for a 13M progenitor in their axisym- 
metric simulations. However, the Oak Ridge group has 
found stron ger and earlier expl osions for a wider range of 
progenitors dBruenn et al. 2009), while in purely Newtonian 
simulations with multi-dimensional neutrino diffusion (in- 
cluding energy dependence but without energy-bin coupling) 
iBurrows et alJ (120061 120071) could not see any success of the 
delayed neutrino-driven mechanism. 

While the reason for the discrepant results of these sim- 
ulations cannot be satisfactorily understood on the basis of 
published results, the marginality of the 2D explosions of the 
Garching group and the lack of neutrino-driven explosions 
in the simulations by IBurrows et alJ (120061 120071) raises the 
important question about the influence of the third spatial 
dimension on the post-bounce evolution of collapsing stel- 
lar cores. Three-dimensional (3D) fluid dynamics with their 
inverse turbulent energy cascade compared to the 2D case 
are likely to change the flow pattern on large scales as well 
as small scales. They could have an influence on the ex- 
istence and the growth rate of nonradial hydrodynamic in- 
stabilities in the super nova core even in th e absence of stel- 
lar rotation (see, e.g.. Ilwakami et all 12008b but in particular 
with a modera te amount of angular momentum in the progen- 
itor s t ar (e.g.. | Blondin & Mezzacap pa 2006a; Iwaka mTet al.l 
120091 iFernandezI 120 101) . and thus could lead to differences 



in the hydrodynamic and thermodynamic conditions for the 
operation of the neutrino-heating mechanism. In particular, 
3D flows might cause important changes of the dwell time of 
postshock matter in the layer where neutrinos deposit energy, 
which is a crucial aspect deciding about the viability and ef- 
ficiency of the neutrino-driven su pernova mechanism (some 
asp ects of this were disc ussed by iMurphv & Burrows! 120081 
an d lMarek & Jankall2009l) . 

Indeed, employing a simplified treatment of neutrino effects 
by including local neutrino-cooling and heating terms for a 
chosen value of the neutrino luminosity and spectral temper- 
ature inst ead of solving the comp utationally intense neutrino 
transport, Nordhaus et al. I (120101) found considerably easier 
and earlier explosions in 3D than in 2D. In the context of 
the concept of a critical value of the neutrino luminosity that 
(for a given mass accretion rate onto the stalled supernova 
shock ) has to be exceeded to obtain neutrino-driven explo- 
sions (IBurrows & Goshvlll993t Uanka & Mullelll996fc Ijankal 
200U lYamasaki & Yamada|l2005t IMurphv & Burrows! l2008h 



Pejcha & Thompsonll201 It lFernandeal2012l) . they quantified 
the improvement of 3D relative to 2D by a 15-25% reduc- 
tion of the critical luminosity value. In particular, they ob- 
served that 3D postshock convection leads to higher average 
entropies in the neutrino-heating layer, thus improving the 
conditions for shock revival due to a significant stretching of 
the residence time of matter i n the layer where it gai ns energy 
from neutrinos. Very recently. [Takiwaki et al.l ( 1201 11) reported 
enhanced maximum dwell times of a small fraction of the ma- 
terial in the gain region in a 3D simulation compared to the 2D 
case of an 1 1 .2 M star, but could not unambiguously link this 
effect to an easier explosion of the 3D model. In particular, 
their 3D simulation showed a shock expansion that was more 
delayed than in the 2D run, and the 3D conditions did not ap- 
pear more favorable for an explosion with respect to a variety 
of quantities like the net heating rate, the heating timescale or 
the profiles of maximum and minimum entropies. 
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In this paper we present a comparative investigation for 
11.2M and 15 M progenitors in one, tw o, and three di- 
mensi ons along the lines of the study by iNordhaus et al. I 
varying the driving neutrino luminosities used in time- 
dependent collapse simulations of the two stars. While our 
results for spherically symmetric (ID) and 2D models basi- 
cally confirm the dimensio n-de pendent differences foun d by 
iMurphv & Burrows! (120081) and INordhaus et"aT1 (120101) , our 
calculations do neither exhibit a strict 1D-2D-3D hierarchy of 
the average entropy in the gain layer, nor do they show any 
clear signs that 3D effects facilitate the development of the 
explosion better than nonradial motions in 2D. Attempting to 
understand the reason for this puzzling finding we vary the 
resolution of the spherical coordinate grid used for our 2D 
and 3D simulations. The outcome of these studies reflects 
the action of the energy flow within the turbulent energy cas- 
cade. The latter transports the driving energy provided by 
neutrino heating and gravitational energy release in the ac- 
cretion flow from small to large scales in 2D and opposite 
in 3D. Models in 2D show growing large-scale asymmetry 
and quasi-periodic time variability and explode clearly easier 
with higher resolution, whereas in 3D better resolved mod- 
els are observed to become more similar to the ID case and 
thus to be farther away from an explosion. This suggests that 
the success of the neutrino-driven mechanism could be tightly 
linked to the initiation of strong non-radial mass motions in 
the neutrino-heated postshock layer on the largest possible 
scales, implying that the easier explosions of our 2D mod- 
els with higher resolution are a consequence of more violent 
activity due to the st anding accretion-shock instability (SASI; 
iBlondin et al]|2003l) . wher eas the better resolved 3D models 
for the employed artificial setup of supernova-core conditions 
tend to reveal considerably reduced amplitudes of low-order 
spherical-harmonics modes of nonradial deformation and thus 
behave more similar to the ID case. 

The paper is structured as follows. In Sect. [2] we briefly 
describe our numerics and implementation of neutrino source 
terms. In Sect. [3] we give an overview of the simulations pre- 
sented in this paper. Our investigations of the dependence of 
the critical luminosity on the spatial dimension will be pre- 
sented in Sect. [4] and results of resolution studies in Sect. [5] 
An interpretation of our findings will follow in Sect. [6] Sec- 
tion [7] contains the summary and conclusions. In App. lAlwe 
present ID simulations that document our efforts to repro- 
duce the results of previous, similar studies in the literature by 
straightforwardly applying the neutrino treatment described in 
these works. 

2. NUMERICAL SETUP 

We solve the equations of hydrodynamics reflecting the 
conservation of mass, momentum, and energy, 



^ + V-(pv) = 0, 
at 

dpy 

+ V • (pv <8> v) + VP = -pVO . 

at 



(1) 

(2) 



de 
8t 



+ V • [(« + P) v] = -pv • V0> +p(Q + v - Q-) , (3) 



where p is the mass density, v the fluid velocity, O the grav- 
itational potential, P the pressure, and e the total (inter- 
nal+kinetic) fluid energy density. These equations are in- 
tegrated in a conservative form (for which reason the en- 



ergy equation is solved for the total energy density) us- 
ing the explicit, finite-volume, higher-or der Eulerian, multi- 
fluid hydrodynamic s code Prometheus dFrvxell et al.l Il99lt 
iMiiiler et alJ 1 1 99 l"allbh . It is a direct i mplementation of the 
Piecew ise Parabolic Method (PPM) of IColella & Woodwardl 
using th e Riemann solver for real gases developed by 
IColella & Glazl d 19851) and the directional splitting approach 
of lStrand (Il968h to treat the multi-dimensional problem. In 
order to prevent odd-even coupling ((Quirk 1994) w e switch 
from t he original PPM method to the HLLE solver o flEinfeldl 
(1 19881) in the vicinity of strong shocks. The advection of nu- 
clear species is treated by the Consistent Mul ti-fluid Advec- 
tion (CMA) scheme of lPlewa & Miillerl d!999t). 

To facilitate comparison with INordhaus' etaiTI d2010h we 
also employ the high-density equation of state (EoS) of 
IShen et all (119981) and do not include general relativistic cor- 
rections. We use the monopole approximation of the Poisson 
equation to treat Newtonian self-gravity. 

To make our extensive para meter study possible, we us e the 
local source terms appl ied bv lMurphv & Burrows! (120081) and 
INordhaus et al. I (120101) instead of detailed neutrino transport 
(see Uankall2001l for a derivation of these source terms). In 
this approach the neutrino heating and cooling rates and 
Qy are given by 



Q + v = 1.544- 10 



2D 



10 52 erg: 



4MeV 



(4) 



100 km\ 2 / \ _ 



erg 



gs 



erg 



(5) 



These approximations depend on local quantities, namely the 
density p, the temperature T, the distance from the center of 
the star r, and the neutron and proton number fractions Y„ and 
Y p , respectively. In Eq. (HJi the electron-neutrino luminosity 
L Vr is a parameter and is assumed to be equal to the electron 
antineutrino luminosity Ly e = L Ve . The neutrino temperature 
T Ye is set to 4 Me V. 

The employed source terms, Eqs. (01 and (0, without the 
factors e~ Tclr are valid for optically thin regions only. In or- 
der to model the transi tion to neutrino t rappin g at high op- 
tical depths, we follow INordhaus et al. I (120101) and multiply 
the heating and cooling terms by e~ TdI to suppress them in the 
inner opaque core of the proto-neutron star. Here, the optical 
depth for electron neutrinos and antineutrinos is defined as 



Teff(r) 



J, ^ 



ff (r')dr'. 



(6) 



The eff ective opacity K e ff wa s derived by iJankal (120011) and 
given in lMurphy et al] (120091) : 



K eff * 1.5-10- 7 -*,, 



10 I0 gcm- 3 



4MeV/ 



cm 



(7) 



where X„ p = \ (y„ + Y p ^j accounts for composition averag- 
ing. In multi-dimensional simulations we evaluate the radial 
integration for the optical depth T^g independently for each 
latitude 6 in 2D and each pair of latitudinal and azimutha l 
angles (8,(f>) in 3D. Note that in IMurphv & Burrows! d2008l) 
the exponential suppression factor is absent in the heating and 
cooling terms (or was not mentioned), which otherwise agree 
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with ours, while no definition of the effective optical depth 
r e ff is given in iNordhaus et al. 1 ( 120101) . The facto r (y„ + Ypj 
in Eq s. (ID and Q is included inlNordhaus et alH (120101). b ut 
not in lMurphv & Burrowsl d2008l) and lMurphv et alj ( 120091) . 

The time period from the onset of the collapse until 1 5 ms 
after bounce was track ed with the Prometheus- Vertex code 
(iRampp & Jankal [20021) in ID including its detailed, multi- 
group neutrino transport, all relevant neutrino reactions and a 
"flashing treatment" for an approximative description of nu- 
clear burning during infall. This means that until 15 ms after 
bounce we describe neutrino effects including the evolution of 
the electron fraction Y e with high sophistication. At 15 ms af- 
ter core bounce we switch to the simple neutrino heating and 
cooling terms and upon mapping from ID to 2D also impose 
(on the whole computational grid) random zone-to-zone seed 
perturbations with an amplitude of 1 % of the density to break 
spherical symmetry. 

Although during the subsequent evolution we apply the 
hea ting and cooling expres sions of Eqs. dH) and (H) follow- 
ing INordhaus eiaLl (120101) and IMurphy & Burrowsl (120081) . 
we refrain from adopting their treatment of c hanges of the 
electro n fraction Y e . Following a suggestion by ILiebendorferl 
(120051) . they prescribed Y e simply as a function of density, 
Y e (p), instead of solving a rate equation with source terms 
for electron neutrino and antineutrino production and destruc- 
tion consistently with the expressions empl oyed for neutrino 
heating and cooling. ILiebendorferl ( 120051) found that such 
a parametrization, supplemented by a corresponding entropy 
source term (and a neutrino pressure term in the equation of 
motion), yields results in good agreement with ID simula- 
tions including neutrino transport during the collapse phase 
until the moment of bounce-shock formation. A universal Y e - 
p-relation, which serves as the basis of this approximation and 
can be inferred for the infalling matter during the homolo- 
gous collapse phase, however, applies neither for the evolu- 
tion of the shocked accretion flow in the post-bounce phase 
nor for expanding neut rino-heated gas (see, e.g., Fig. 4.9 in 
iThielemann et al.11201 11) . For example, a comparison with su- 
pernova models with detailed neutrino transport shows that 
the F e -p-relation fitted to the homologous phase overestimates 
the deleptonization of the gas in most of the gain layer but un- 
derestimates the lepton loss of matter in the cooling layer and 
neutrinospheric region. Moreover, the question arises how the 
lepton evolution shall be treated in matter that reexpands and 
thus moves from high to low density? Even more, the entro py 
source term introduced in Eq. (5) of ILiebendorferl (120051) is 
designed to specifically account for gas-entropy changes due 
to neutrino production by electron captures and subsequent 
energy transfers in (neutrino-electron) scatterings. The cor- 
responding energy loss or gain rate of the medium through 
the escape or capture of electron neutrinos with mean energy 
E Ve , which is given by 5Q v JSt - E Ve 6Y e /6t, is not i nclud ed 
in the heating and cooling terms Q+ and Q~ of Eqs. (1415 \ of 
the present work. Adding it as an extra term would imply par- 
tial double-counting of the energy exchange via electron neu- 
trinos, and omitting it means to overestimate the entropy in- 
crease in infalling, deleptonizing matter and to underestimate 
entropy gains of decompressed gas with growing Y e . Because 
of the long list of such inconsistencies, whose implications 
are hard to judge or control, we decided to ignore Y e changes 
of the stellar medium in our post-bounce simulations (except 
for the models discussed in App.lAl. 

This choice can be justified, but it is certainly not a perfect 



approach because it may also exclude effects of importance 
in real supernova cores, whose physical processes require the 
inclusion of neutrino transport for an accurate description of 
the energy and lepton-number evolution. One of the unde- 
sired consequences of keeping Y e fixed in the accretion flow 
after core bounce is an overestimation of the electron pressure 
in the gas settling onto the forming neutron star. In order to 
enforce more compression and thus to ensu re close similar- 
ity of o ur re sults to the ID mode l s stud ied bv lNordhaus et al.l 
(120101) and IMurphy & Burrowsl (120081) . we have to enhance 
the net cooling of the accreted matter by reducing the effec- 
ti ve opacity K e g by a fa ctor of 2.7 compared to the value given 
in lMurphv et al.l (120091) and in Eq. (7). This reduction factor is 
chosen such that our simulations reproduce the minimum val- 
ues of the critical luminosity found to be necessary for trigger- 
ing explosions of the 15 M progenitor (2.6 ■ 10 52 erg s _1 ) and 
fo r the 11.2 M G star ( 1 . 3 • 1 52 erg s~ 1 ) in the ID simulations 
o flMurphv & Burrowsl ( 120081) . Without the reduction factor of 
K e ff , our models turn out to explode too easily because of weak 
cooling (see the results in App. [A}. We stress that any expo- 
nenti al su ppression factor of the heating and cooling rates in 
Eqs. ( 1415b is a pragmatic and ad hoc procedure to bridge the 
transition from the optically thin to the optically thick regime, 
where neutrino transport is most complicated. From transport 
theory neither the exponential factor nor the exact definition 
of the optical depth of the exponent can be rigorously derived. 

In our reference set of standard simulations, we employ 400 
non-equidistant radial zones, which are distributed from the 
center to an outer boundary at 9000 km. The latter is suf- 
ficiently far out to ensure that the gas there remains at rest 
during the simulated evolution periods. The radial zones are 
chosen such that the resolution Ar/r is typically better than 
2%. In the multi-dimensional models with standard resolu- 
tion, we employ a polar grid with an angular bin size of 3° (60 
9- and 120 0-zones). For the high-resolution ID and multi- 
dimensional models of Sects. |4]and|5]we compute with up to 
800 radial zones and in 2D with an angular resolution down to 
0.5°, in 3D down to 1.5°. The additional radial grid zones are 
distributed such that the region between 20 and 400 km, i.e., 
not only the cooling layer around the neutrinosphere but also 
the gain layer between gain radius and shock, are resolved 
significantly better. While full convergence of the ID results 
requires 600 radial zones or more, most of the plots compar- 
ing models for different dimensions show simulations with 
our standard resolution of only 400 radial cells (unless stated 
otherwise), because this is the limit for which we could per- 
form a larger set of 3D runs in acceptable time. To avoid an 
extremely restrictive Courant-Friedrich-Lewy (CFL) timestep 
in our 3D calculations, we simulate the inner core above a 
density of p — 10 12 g/cm 3 in spherical symmetry. 

3. INVESTIGATED PROGENITORS AND MODELS 

Our models are base d on th e 15 M progenitor star sl5s7b2 
of iWooslev & Weaver! d 19951) and an 11.2 M Q progenitor of 
IWooslev et al.l (120021) . The calculations for these progenitors 
with our standard angular resolution of 3° are summarized in 
TableQ] This table is arranged such that horizontal rows have 
the same driving luminosities for simulations performed in 
different dimensions. Varying the prescribed driving luminos- 
ity L Ye from run to run we present for each of the 11.2 M and 
15 M progenitors several ID, 2D, and 3D simulations. All 
of the ID and 2D simulations cover at least Is after bounce. 
The nonexploding 3D simulations with standard angular res- 
olution of 3 degrees were not stopped until at least 600 ms af- 
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Fig. 1. — Time evolution of the mass accretion rate, M(r) = Anr 1 p(r)\v(r)\, 
evaluated at 500 km for the 1 1.2 M Q and the 15 M Q progenitors in nonexplod- 
ing models. 



TABLE 1 

1 1.2 Mq AND 15 M Q RESULTS WITH STANDARD GRID OF 400 RADIAL 
ZONES AND 3 DEGREES ANGULAR RESOLUTION. 



ID 2D 3D 

J a t b M c t M t M 

M* e i exp J ' J exp <exp iw exp 'exp J "exp 

(10 52 erg/s) (ms) (M /s) (ms) (M G /s) (ms) (M G /s) 



11.2 Mr 



0.8 














0.9 










731 


0.085 


1.0 






563 


0.082 


537 


0.086 


1.1 






461 


0.091 






1.2 






357 


0.104 


319 


0.112 


1.3 


819 


0.082 


307 


0.114 






1.4 


499 


0.088 


241 


0.126 


232 


0.130 


1.5 


380 


0.100 


232 


0.130 






1.6 


345 


0.106 


203 


0.137 






15 M e 


2.0 














2.1 














2.2 






876 


0.197 


612 


0.226 


2.3 






428 


0.261 


426 


0.261 


2.4 






442 


0.255 






2.5 






283 


0.313 


281 


0.314 


2.6 


710 


0.215 


285 


0.312 






2.7 


489 


0.247 


262 


0.316 






2.8 


390 


0.271 


242 


0.322 


236 


0.324 


2.9 


281 


0.314 


235 


0.325 






3.0 


258 


0.318 


236 


0.324 






3.1 


248 


0.320 


220 


0.327 







d Electron-neutrino luminosity. 



Time after bounce of onset of explosion. A "— " symbol indicates that the 
model does not explode during the simulated period of evolution. 
c Mass accretion rate at onset of explosion. 

ter bounce. For models with higher resolution the simulation 
times are given in Table [2] 

In TableQ]the time of the onset of an explosion, f exp , and the 
mass accretion rate at that time, M e xp, are listed as characteris- 
tic quantities of the models. The beginning of the explosion is 
defined as the moment of time f exp when the shock reaches 
an average radius of 400 km (and does not return lateron), 
while nonexploding models are denoted by a "-" symbol. In 
multidimensional simulations the corresponding shock posi- 
tion is defined as the surface average over all angular direc- 
tions, (Rs) = §dQRs(Q,). The lowest driving luminosity 



yielding an explosion for a given va lue of the mass accretion 
rate is termed the critical lu minosity dBurrows & Goshvl 19931: 
iMurphv & Burro ws 2008). We determine the mass accretion 
rate M(r) = 47ir 2 p(r) |v(r)| at the time of the onset of the ex- 
plosion just exterior to the shock, i.e. at a radius of 500 km, 
where the infalling envelope is spherical (except for the small 
seed perturbations imposed on the multi-dimensional mod- 
els). In Fig. Q]the mass accretion rates of the 11.2M and 
15M Q progenitors are depicted for the nonexploding runs 
with the lowest driving neutrino luminosities. Because the 
shock can be largely deformed in multi-dimensional simula- 
tions and its outermost parts can extend beyond 500 km (and 
thus impede a clean determination of the mass accretion rate) 
when its average radius just begins to exceed 400 km, we refer 
to the functions plotted in Fig.[T]for defining the mass accre- 
tion rate at the time when the explosion sets in. 

4. CRITICAL LUMINOSITY AS FUNCTION OF DIMENSION 

Based on steady-state solutions of neutrino-heated and 
-cooled accretion flows be tween the stalled s hock and the 
proto-neutron star surface, iBurrows & Goshvl d 1993b identi- 
fied a critical condition that can be considered to separate ex- 
ploding from nonexploding models. They found that for a 
given value of the mass infall rate M onto the accretion shock 
steady-state solutions cannot exist when the neutrino-heating 
rate in the gain layer is sufficiently large, i.e., for neutrino 
luminosities above a threshold value L v . This result can be 
coined in terms of a critical condition L y (M) expressing the 
fact that either the driving luminosity has to be sufficiently 
high or the damping mass accretion rate enough low for an 
explosion to b ecome possible. The crit ical M-L v -curve was 
interpreted by IBurrows & Goshvl fi993j) as a separating line 
between the region above, where due to the non-existence 
of steady-state accretion solutions explosions are expected to 
take place, from the region below, where neutrino energy in- 
put behind the shock is insufficient to accelerate the stalled 
shock outwards and thus to trigger an explosion. 

This interpretation of the steady-state results was consis- 
tent with hydrodynamical sim ulations of collapsing a nd ex- 
ploding stars in ID and 2D by I Janka & Mullerl (119961) . Per- 
forming a more extended p arameter study than the latter work, 
IMurphv & Burrows! (120081) explored the concept of a critical 
condition systematically with time-dependent hydrodynami- 
cal models. They showed that a critical luminosity indeed 
separates explosion from accretion and confirmed that this 
value is lowered by ~30% when going from spherical sym- 
metry to two dimensions, at least when a fixed driving lumi- 
nosity is adopted and feedback effects of the hydrodynamics 
on the neutrino emission are ignored. Some 2D effects includ- 
ing possible consequences of rotation were discussed before 
on grounds o f steady -state models by lYamasaki & Yamadal 
(I2005L 120061) . while Uankal (1200 lh tried to include time- 
dependent aspects of the shock-revival problem and took 
into account an accretion component of the neutrino lumi- 
nosity in addition to the fixed core component. The influ- 
ence of such an accretion contribution was more recently 
also estimated by iPeicha & Thompson! d201 ll) . who solved 
the one-dimensional steady-state accretion problem between 
the neutron star a n d the accretion shock along the lines of 
IBurrows & Goshvl (119931) . but attempted to obtain a deeper 
understanding of the critical condition by comprehensively 
analysing the structure of the accretion layer and of the lim- 
iting steady-state solution in dependence of the stellar condi- 
tions. They found that the critical value for the neutrino lumi- 
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Fig. 2. — Critical curves for the electron-neutrino luminosity (L Ve ) versus mass accretion rate (M) (left plot) and versus explosion time i cxp (right plot) for 
simulations in ID (black), 2D (blue), and 3D (red) with standard resolution. The accretion rate is measured just outside of the shock at the time ? cxp when the 
explosion sets in. For the 15 M Q progenitor ID results are displayed for 400, 600, and 800 radial zones. Higher radial resolution compared to the standard ID 
runs makes explosions slightly more difficult; convergence is achieved for >600 radial bins. Multi-dimensional results are shown for different angular resolutions, 
where available, but always computed with 400 radial zones. Note that 2D simulations with improved angular zoning explode more easily, whereas in 3D only 
one case was computed (the 11.2 M Q simulation for L Ve = 1 .0 ■ 10 52 erg s _I ) that developed an explosion also with better angular resolution. 
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Fig. 3 . — Time evolution of the average shock radius as function of the post- 
bounce time /pb for simulations in one (thin dashed lines), two (thin solid 
lines), and three dimensions (thick lines). The shock position is defined as 
the surface average over all angular directions. The top panel shows results 
for the 11.2 M B progenitor and the bottom panel for the 1 5 M Q progenitor, all 
obtained with our standard resolution. Different electron-neutrino luminosi- 
ties (labelled in the plots in units of 10 52 erg s~') are displayed by different 
colors. 

nosity is linked to an "antesonic condition" in which the ratio 
of the adiabatic sound speed to the local escape velocity in the 
postshock layer reaches a critical value above which steady- 
state solutions of neutrino-heated accretion flows cannot be 
obtained. By performing high-resolution hydrodynamic sim- 



ul ations [Fernand ez ( 20 1 2) found that radial instability is a suf- 
ficient condition for runaway expansion of an initially stalled 
core-collapse supernova shock if the neutrinospheric param- 
eters do not vary with time and if heating by the accretion 
luminosity is neglected. However, the threshold neutrino lu- 
minosities for the transition to runaway instability are in gen- 
eral different from the limitin g values for steady - state solu- 
tions of the kind discussed bv lBurrows & Goshvl d 1 993b and 
iPejcha & Thompson! d201 lb . iNordhaus et al. I (120101) general- 
ized th e hydrodynamic investigations of iMurphy & Burrows! 
(120081) to include 3D models and found another reduction of 
the threshold luminosity for explosion by 15-25% compared 
to the 2D case. 

Despite the basic agreement of the outcome of these investi- 
gations it should be kept in mind that it is not ultimately clear 
whether the simple concept of a critical threshold condition 
separating explosions from failures (and the dependences of 
this threshold on dimension and rotation for example) holds 
beyond the highly idealized setups considered in the men- 
tioned works. None of the mentioned systematic studies by 
steady-state or hydrodynamic models was able to include ad- 
equately the complexity of the feedback between hydrody- 
namics and neutrino transport physics. In particular, none of 
these studies could yield the proof that the non-existence of 
a steady-state accretion solution for a given combination of 
mass accretion rate and neutrino luminosity is equivalent to 
the onset of an explosion. The latter requires the persistence 
of sufficiently strong energy input by neutrino heating for a 
suffientl y long period of time Thi s is especially important 
because lPejcha & Thompso"n](l201 lb showed that the total en- 
ergy in the gain layer is still negative even in the case of the 
limiting accretion solution that corresponds to the critical lu- 
minositjO]. Within the framework of simplified modeling se- 
tups, however, the question cannot be answered whether such 
a persistent energy input can be maintained in the environ- 
ment of the supernova core. 

Following the previous investigation s by 

IMurphy & Burrows! d2008h and INordhaus et aTl d201(X) 
we performed hydrodynamical simulations that track the 



1 Note that Fernandez 12012) demonstrated that transition to runaway oc- 
curs when the fluid in the gain region reaches positive specific energy. 



6 



Hanke, Marek, Miiller, & Janka 



post-bounce evolution of collapsing stars for different, fixed 
values of the driving neutrino luminosity. Since the mass 
accretion rate decreases with time according to the density 
profile that is characteristic of the initial structure of the 
progenitor core (see Fig. Q] for the 11.2 and 15 M Q stars 
considered in this work), each model run probes the critical 
value of M exp at which the explosion becomes possible for 
the chosen value of L v = L Ve = Ly e . The collection of 
value pairs (M exp ,L Vf ) defines a critical curve L V {M). These 
are shown for our ID, 2D, and 3D studies with standard 
resolution for both progenitor stars in the left panel of Fig. [2] 
and in the case of the 15 M st a r can be directly compared 
with Fig. 1 of iNordhaus et al71 (120101) . Table [J lists, as a 
function of the chosen L Ve , the corresponding times f exp when 
the onset of the explosion takes place and the mass accretion 
rate has the value of M exp . The post-bounce evolution of a 
collapsing star proceeds from high to low mass accretion rate 
(Fig. Q3, i.e., from right to left on the horizontal axis of the 
left panel of Fig. [2] When M reaches the critical value for 
the given L v , the model develops an explosion. The right 
panel of Fig.f2]visualizes the functional relations between the 
neutrino luminosities L Ve and the explosion times f e xp for both 
progenitors and for the simulations with different dimensions. 

At first glance Fig. [2] r eproduces bas i c trend s that are vis- 
ib le in Figs. 17 of IMurphv & Burrowsl (12008b and in Fig. 1 
of INordhaus et al. 1 (120101) . For example, the critical luminos- 
ity increases for higher mass accretion rate and the values for 
spherically symmetric models are clearly higher than those 
for 2D simulations. However, a closer inspection reveals in- 
teresting differences compared to the previous works. 

• In general the slopes of our critical L Vf (M)-relations 
appear to be considerably steeper and in the case of 
the 15M Q star they exhibit a very steep rise above 
M =s O.31M s _1 . This means that explosions for 
higher mass accretion rates are much harder to obtain 
and therefore the tested driving luminosities in our sim- 
ulations do not lead to explosions earlier than about 
200 ms after bounce, independent of whether the mod- 
eling is performed in ID, 2D or 3D. In cont rast, the 
critical curves given by INordhaus et al.l (120101) show a 
moderately steep increase over the whole range of plot- 
ted mass accretion rates between about 0.1 M s^ 1 and 
more than 0.5 M Q s . 

• Nonradial flows in the 2D case, by which the residence 
time of accreted matter in the gain layer could be ex- 
tended or more matter could be kept in the gain re- 
gion, reduce the critical luminosities by at most ~15% 
of the ID- values for the 15M star and <25% for 
the 11.2 M model, w hich is a somewhat s maller dif- 
ference than found by IMurphv & Burro ws (2008)) and 
INordhaus et al. I dMIO) . 

• Most important, howev er, is the fact that we can not con- 
firm the observation bv lNordhaus et al. I (120101) that 3D 
provides considerably more favorable conditions for an 
explosion than 2D. Our critical curves for the 2D and 
3D cases nearly lie on top of each other. There are 
minor improvements of the explosion conditions in 3D 
visible in both panels of Fig. [2] and the numbers of Ta- 
ble[T] e.g., a 10% reduction of the smallest value of L v 
for which an explosion can be obtained for the 11.2 M 
star, a ~ 260 ms earlier explosion for the lowest lumi- 
nosity driving the 15 M Q explosion (2.2 ■ 10 52 ergs~'), 
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Fig. 4. — Time evolution of the mass-weighted average entropy in the gain 
region for one-dimensional (thin dotted lines), two-dimensional (thin solid 
lines), and three-dimensional (thick lines) simulations with different angular 
resolutions (corresponding to different colors). The top panel displays the 
11.2 M Q results for an electron-neutrino luminosity of L Ve = 1 .0- 10 52 erg s , 
the middle panel shows the 15 M Q runs for an electron-neutrino luminosity 
of Ly = 2.1 ■ 10 52 erg s , and the bottom panel the 15 M Q models for Ly e = 
2.2 ■ 10 52 erg s . The strong decrease of the average entropy that terminates 
a phase of continuous, slow increase signals the onset of the explosion when 
a growing mass of cooler (low-entropy) gas is added into the gain layer after 
being swept up by the expanding and accelerating shock wave. 

and a tendency of slightly faster 3D explosions for all 
tested luminosities (see Fig.|3]and Table[T]i. All of these 
more optimistic 3D features, however, will disappear 
for simulations with higher resolution as we will see in 
Sect. [5] 

Before we discuss the o rigin of the differen c es bet ween 
our results and those of IMurphv & Burrowsl (120081) and 
INordhaus et al. I (120101) we would like to remark that the kinks 
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Fig. 5. — Scatter-plots of the entropy structure as function of radius for simulations of the 11.2M progenitor with an electron-neutrino luminosity of L Ve = 
1.0- 10 52 ergs-' at 400 ms (left) and 600 ms (right) after core bounce. The red dots correspond to the 2D results, black ones to 3D, the light blue line is the 
entropy profile of the ID simulation, and the dark-blue and green curves are mass-weighted angular averages of the 2D and 3D models, respectively. Both 
multi-dimensional simulations were performed with an angular resolution of two degrees and both yield explosions (at ~530ms p.b. in 2D and ~570ms p.b. in 
3D; see TableffJ- Note that different from Fig. [3] unshocked material at a given radius is included when computing angular averages. The dispersion of entropy 
values in the uns hocked flow of 2D a nd 3D simulations is a consequence of the imposed density-seed perturbations (cf. Sect.ff), which grow in the supersonical 
infall regime (see Buras et al. 2006b). 
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Fig. 6. — Scatter-plots of the entropy structure as function of radius for simulations of the 15M Q progenitor with an electron-neutrino luminosity of L Ve = 
2.1 • 10 52 ergs-' at 350 ms (left) and 700 ms (right) after core bounce. The red dots correspond to the 2D results, black ones to 3D, the light blue line is the 
entropy profile of the ID simulation, and the dark-blue and green curves are mass-weighted angular averages of the 2D and 3D models, respectively. Both 
multi-dimensional simulations were performed with an angular resolution of 1.5 degrees. While the 2D model develops an explosion setting in ~720ms after 
bounce, the 3D model does not produce an explosion (TableffJ- Note that different from Fig. [4] unshocked material at a given radius is included when computing 
angular averages. The dispersion of entropy values in the uns hocked flow of 2D a nd 3D simulations is a consequence of the imposed density-seed perturbations 
(cf. Sect. |2), which grow in the supersonical infall regime (see Buras et al. 2006b). 
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and even nonmonotonic parts of the curves shown in Fig. [2] 
in particular for the multi-dimensional cases are connected 
to our definition of the explosion time as being the moment 
when the mean shock radius exceeds 400 km. Especially in 
cases where the shock deformation is large (which is an issue 
mainly in some of the 2D simulations) this definition is asso- 
ciated with significant uncertainty in the determination of the 
exact explosion time f exp and therefore also of the correspond- 
ing mass accretion rate M exp . 

In addition to the results for our standard resolution, Fig. [2] 
presents ID models with higher radial resolution (600 and 
800 radial zones). The critical curves with better zoning ex- 
hibit the same overall slopes as those of the standard runs, 
but there is a slight shift towards higher values of the critical 
luminosity (or, equivalently, a small shift to lower M exp and 
later f exp ). This is caused by somewhat larger neutrino energy 
losses from the cooling layer with better radial resolution, an 
effect which makes explosions more difficult. We stress that 
this resolution-dependent cooling effect is a consequence of 
the simplified neutrino-loss treatment. The employed cool- 
ing rate is not able to reproduce real transport behavior and 
leads to the development of a pronounced gaussian-like den- 
sity peak (with density excess of a few compared to its sur- 
roundings) in the neutrino-cooling layer. This local density 
maximum has a strong influence on the integrated energy loss 
by neutrino emission. With better zoning the peak becomes 
better resolved and even grows in size. Convergence of ID 
results seems to be achieved for >600 radial zones, but in 2D 
and 3D the artificial density peak prevents numerical conver- 
gence for all employed radial zonings (cf. Sect. [5). Figure [2] 
also displays some data points for multi-dimensional models 
with better angular resolution (all of them, however, com- 
puted with 400 radial mesh points). These will be discussed 
in Sect. [5] 

A more detailed analysis, which we will report on below 
and in App.lAl reveals that the exact values of the critical lu- 
minosities as well as the detailed slope of the critical curves 
seem to depend strongly on the employed description of neu- 
trino effects, whose implementation is subject to a significant 
degree of arbitrariness if detailed neutrino transport is not in- 
cluded i n the modeling (cf. the dis cussion in Sect. |2). The 
fact that iMurphv & Burrowsl (120081) found fairly good over- 
all agreeme nt between their cri t ical re lations L V {M) and those 
obtained by lBurrows & Goshyl (119931) is likely to be linked to 
a basically similar treatment of the neutrino effects. 

Before ~0.2 s after bounce the mass accretion rate in the 
case of the 15 M progenitor changes much more rapidly than 
during the subsequent evolution (Fig. UJ. For the correspond- 
ing M values in excess of ~0.31 M s _1 , the accretion shock 
is therefore not as close to steady-state conditions as later on. 
We see a steep rise of our critical curves at M > 0.31 M s _I 
(fexp 55 0.25 s), which i s a very prominen t differ ence com- 
pared to the results of iBurrows & Goshvl (119931) . who as- 
sumed steady-state accretion, bu t in particular also compare d 
to t he hydrodynamic resul ts of IMurphv & Burrowsl (120081) . 
and iNordhaus et afl (1201 Oh even in the ID case. In order to 
explore possible reasons for this difference, we performed ID 
simulations in which t he neutrino treatment is copied from 
INordhaus et al. I (120101) as closely as possible (i.e., the re- 
duction factor of 2.7 in the exponent of e~ T °" is not applied 
and deleptonization is taken into account by using a Y e (p) re- 
lat ion, but not the cor responding entropy changes proposed 
bv lLiebendorfenl2005l: for more details on these results, see 



App. [Aj. These runs reveal that the steep rise of our L V (M)- 
curves is caused by a strong increase of the neutrino-cooling 
rate with higher values of M, in particular when we apply 
our neutrino treatment. The corresponding energy losses in- 
hibit explosions for low values of the driving luminosity. The 
stronger cooling is linked mainly to our reduction of the ef- 
fective optical depth T e ff, which we had to apply in order to 
reconcile the mass accretion rates and e xplosion times with 
the low est driving luminosities for which lMurphv & Burrowsl 
( 120081) had obtained explosions for the 1 1.2 and 15 M Q stars 
(cf. Sect. |2j. For example, in the case of the 15 M progen- 
itor a driving luminosity of L Vg - 3.1 ■ 10 52 ergs~' triggers 
an explosion at fexp ~ 250ms p.b. and M exp * 0.32 M s~' 
(Table[T]and Fig. [2J|, where as with an implementat ion of neu- 
trino effects closer to that of Nordh aus et al. 1(120101) the explo- 
sion sets in at f™ =s 120 ms p.b. and M exp ~ 0.8M o s _1 (see 
Fig.[19]in App.lAl. Shortly before this moment (at 75 ms after 
bounce) the total energy loss by neutrino cooling is abo ut 10 
times lower with the scheme of INordhaus et al. T (12010b than 
with our neutrino implementation. The latter yields an inte- 
grated energy-loss rate of ~9T0 52 erg s _1 and significant cool- 
ing even at dens i ties be tween 10 12 and 10 13 g cm -3 , where the 
INordhaus et al. I (120101) treatment shows essentially no cool- 
ing. Neither the magnitude of the total neutrino-energy loss 
rate nor the region of energy extraction with our modeling 
approach are implausible and in disagreement with detailed 
transport simulations during a stage when the ma ss accretion 
rate st ill exceeds 1 M H s' 1 (cf., e.g., Fig. 20 in iBuras et al.l 
I2006al) . In contrast, the INordhaus et alH (120101) treatment ap- 
pears to massively underestimate the neutrino energy extrac- 
tion from the accretion flow during this time. 

These findings demonstrate that the results of the critical 
L Vf (M)-relation in ID can be quantitatively as well as qualita- 
tively different with different approximations of neutrino heat- 
ing and in particular of neutrino cooling. Moreover, this gives 
reason for concern that the differences of the critical explosion 
condit ions for 2D and 3D simulations seen bv lNord haus et al. 
( 120101) might have been connected to their treatment of the 
neutrino physics, in particular also because the decrease of 
the critical luminosity from 2D to 3D they found was only 
15-25%, which is a relatively modest change (much smaller 
than the 1D-2D difference) and thus could easily be overruled 
by other effects. Our results for 2D and 3D simulations with 
a different implementation of neutrino sources confirm this 
concern. 

In the INordhaus et aTl (120101) paper the average entropy of 
the matter in the gain region, (s(f)), was considered to be a 
suitable diagnostic quantity that reflects the crucial differenes 
of ID, 2D, and 3D simulations concerning their relevance for 
the supernova dynamics. While in the spherically symmet- 
ric case accreted matter moves through the gain region on the 
shortest, radial paths, nonradial motions can increase the time 
that shock-accreted plasma can stay in the gain layer and ab- 
sorb energy from neutrinos. This can raise the mean entropy, 
internal energy, and pressure in the postshock region and thus 
support the revival of the stalled supernova shock . This seems 
to happen in the simulations bv lNordhaus et al.l (120101) . who 
found that turbulent mass motions in 3D can even improve the 
conditions for the neutrino-heating mechanism compared to 
the 2D case. A crude interpretation of this difference can be 
given by means of random-walk arguments, considering the 
mass motions in convective and turbulent structures as diffu- 
sive process in the gain layer dMurphv & Burrowsl2008l) . The 



SASI activity as key to successful neutrino-driven SN explosions? 



9 



question, however, is whether this effect is a robust 2D-3D 
difference and whether it is the crucial key to successful ex- 
plosions by the neutrino-heating mechanism. 

Our results at least raise doubts. Figure|4]displays the time 
evolution of the mean entropy in the gain layer for 1 1 .2 and 
a 15 M Q models computed with driving luminosities near the 
minimum value for which we obtained explosions. While the 
1 1.2 M Q model explodes with a luminosity of 1 .0 ■ 10 52 erg s 
for all tested resolutions in more than one dimension despite 
minimal differences between the values of {s(t)) compared to 
the ID counterpart, the 15 M progenitor develops an explo- 
sion for the chosen luminosity of 2.1 ■ 10 52 erg s _1 only in the 
case of higher-resolution 2D runs (this will be further dis- 
cussed in Sect. [5]). These successful cases, however, do not 
stick out by especially high values of (s(t)). On the contrary, 
they even have lower mean entropies than the unsuccessful 
3D models! It is obvious that Fig.|4]does not exhibit the clear 
1D-2D -3D hierarchy visible in Fig. 5 of the [Nordhau s et alTI 
(l2OT0l) paper, which was found there to closely correlate with 
the explosion behavior of their models. Instead, the differ- 
ences between simulations in the different dimensions are 
fairly small, and even two-dimensional flows, which unques- 
tionably allow for explosions also when none happen in ID, 
do not appear more promising than the ID case in terms of 
the average entropy of the matter in the gain layefl Similarly, 
3D models possess slightly (insignificantly?) higher values of 
(s(t)) but do not show a clear tendency of easier explosions, 
in particular not the better resolved simulations (see Sect.[5]l. 

The radial entropy structures of ID, 2D, and 3D runs for 
both progenitors, shown in Figs.[5]and[6]once before an explo- 
sion begins and another time around the onset of an explosion 
in at least one of the runs, demonstrate that low-entropy accre- 
tion downdrafts and high-entropy rising plumes of neutrino- 
heated plasma lead to large local variations of the entropy per 
baryon of the matter in the gain layer (scatter regions in the 
plots). However, the mass-weighted angular averages of the 
entropies reveal much smaller differe nces between the I D and 
2D cases than visible in Fig. 4 of thelNordhaus et al. I ( 120101) 
paper and in Fig. 13 of [Murphy & Burrows! (120081) . and ex- 
hibit no obvious signs of more advantageous explosion con- 
ditions in the 3D cases compared to 2D. The noticeable dif- 
ferences in the radial profiles seem to be insufficient to cause 
major differences in the mean entropies computed by addi- 
tional radial averaging (see Fig.|4|. 

How ca n this discrepancy compa red to iNordhaus et al. I 
d2010l) and lMurphv & Burrowsl (120081) be explained, and how 
can one understand the fact that 2D effects play a supportive 
role for neutrino-driven explosions? We will return to these 
questions in Sect. [6j but before that we shall present our re- 
sults of multi-dimensional simulations with varied resolution 
in the following section. 

5. MODELS WITH HIGHER RESOLUTION 

In order to test whether our results for the multi- 
dimensional models depend on the agreeably moderate 3° an- 
gular resolution used in the standard runs, we performed a 
large set of simulations with finer grid spacing especially in 
the angular directions, but also in radial direction. For this 

2 We stress that our basic findings are independent of the exact way how the 
gain radius of the multi-dimensional models is determined, i.e., whether the 
evaluation is performed with an angularly averaged gain radius or a direction- 
dependent gain radius. The outer boundary of the integration volume is de- 
fined by the shock position, which usually forms a non-spherical surface in 
the multi-dimensional case. 



TABLE 2 

Multidimensional models with different resolution. 
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Fig. 7. — Evolution of the average shock radius as a function of time (in sec- 
onds after bounce) for one-dimensional (dashed lines), two-dimensional (thin 
solid lines), and three-dimensional (thick solid lines) simulations employ- 
ing different angular resolutions (color coding). The top panel displays the 
11.2 M Q model for an electron-neutrino luminosity of L,^ = 1 . • 1 52 erg s~ 1 , 
the middle panel shows the 15 M Q star for an electron-neutrino luminos- 
ity of L Ve = 2.1 ■ 10 52 ergs~', and the bottom panel the 15 M Q results for 
L Ve = 2.2- 10 52 ergs-'. 

purpose we concentrated on cases around the minimum val- 
ues of the driving luminosity that triggered explosions of both 
progenitors in our standard runs. The results are listed in Ta- 
ble [2] They indicate a very interesting trend: 2D models with 
finer angular zoning tend to explode more readily, whereas 
better angular resolution in 3D simulations turns out to have 
the opposite effect. 

In the case of the 11.2M progenitor, for example, the 3D 
explosion found to set in at about 730 ms p.b. with an angular 
zone size of 3° for L Ve = 0.9 • 10 52 ergs~' cannot be repro- 
duced with an angle binning of 2°. Moreover, a luminosity of 
L Ve = 1.0 ■ 10 52 erg s _1 leads to an explosion of the 3D model 




Fig. 8. — Evolution of the standard deviation for the shock asphericity as 
a function of post-bounce time for two-dimensional (thin solid lines) and 
three-dimensional (thick solid lines) simulations employing different angular 
resolutions (color coding). As in Fig.UJthe top panel displays the exploding 
1 1.2 M Q models for an electron-neutrino luminosity of L Vl , = 1.0- 10 52 erg s _I , 
The middle panel contains the results for the 15M Q star with an electron- 
neutrino luminosity of L Ve = 2.1 ■ 10 52 ergs -1 , where 2D runs with higher 
resolution lead to explosions while 3D runs do not. The bottom panel shows 
the 15 M Q case for L Ve = 2.2 ■ 10 52 erg s _I . It is remarkable that the 3D run 
with 2° angular resolution does not explode whereas the one with angle bins 
of 3° explodes earlier than its 2D counterpart and develops a very large shock 
deformation at the time the explosion sets in. 

at ~540 ms after bounce with a 3°-grid, but ~35 ms later when 
2° are used. The corresponding 2D models show the inverse 
trend as visible in the top panel of Fig. [7] 

Note that the average shock radii plotted in Fig. [7] as well 
as Fig. [3] exhibit alternating periods of increase and decrease 
in particular in 2D simulations. Such features are a conse- 
quence of the strong sloshing motions of the accretion shock 
and of the associated time-dependent, large global shock de- 
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Fig. 9. — Snapshots of the evolution of the 11.2 M a model with an electron-neutrino luminosity of L Ve = 1.0 ■ 10 52 ergs -1 and 2° angular resolution at post- 
bounce times of ? p b = 300, 350, 400, and 600 ms (from top left to bottom right). The color coding represents the entropy per nucleon of the stellar plasma. The left 
half of each panel displays the entropy distribution for a 2D (axisymmetric) simulation, the right half shows the structure in a meridional cut of the corresponding 
3D simulation. Both models explode after roughly 550 ms after bounce (see Fig. [7] top panel, and Table[2j. Note that the structures of low-entropy downdrafts 
and high-entropy plumes in the neutrino-heating region are rather similar for both cases, despite considerably larger SASI sloshing motions of the shock and 
postshock layer in 2D. When the explosion has set in, the 2D model exhibits an apparent prolate deformation whose development is supported by the symmetry 
axis defining a preferred direction of the 2D system. While the 3D explosion does not appear to be as strongly distorted (in particular the shock surface looks 
more spherical), the post shoc k flow in this case also develops a pronounced hemispheric (dipolar) asymmetry, which can be more clearly seen in the upper and 
lower right panels of Fig. |10| 



formations, which are typical of violent activity by the SASI. 
In 3D the corresponding wiggles and local maxima of the av- 
erage shock trajectory are much less pronounced. A measure 
of the degree of shock asphericity, irrespective of the rela- 
tive weights of different spherical harmonics components, is 
the standard deviation of the shock radius defined by <xs = 



4 



^dD.[Rs(D) - (Rs)] 2 - The standard deviations corre- 
sponding to the average shock radii of Fig. [7J are plotted in 
Fig. [8] which confirms the mentioned difference between 2D 
and 3D runs. 

In spite of this 2D-3D difference of the shock aspheric- 
ity, an inspection of cross-sectional snapshots of our best re- 
solved simulations of the 1 1.2M progenitor with an electron- 
neutrino luminosity of L Vr = 1.0 ■ 10 52 ergs _1 reveals that 
the sizes of the convective plumes and the structure of the 
neutrino-heated postshock layer are fairly similar in the 2D 
and 3D cases before explosion (which in both models de- 
velops shortly after 500 ms): In Fig. [9] it is difficult to judge 
by eye inspection whether the displayed simulation was con- 
ducted in 2D (left half -panels) or 3D (right half-panels). Even 



after the explosions have taken off the global deformation of 
the shock in both cases is not fundamentally different in the 
sense that low-order spherical harmonics modes (dipolar and 
quadrupolar components) determine the global asymmetry of 
the shock surface and in particular of the distribution of down- 
drafts and expanding bubbles in the gain region (see Fig. [9] 
lower right panel, and Fig.[T0il. At a closer inspection one can 
notice some secondary differences in the morphology of the 
convective and downflow features. Despite the same angu- 
lar resolution the images of Figs. l9l and [TOl reveal more small 
structures in the 3D case compared to the 2D data, which ap- 
pear more coherent, smoother, and less fragmented into finer 
substructures and filaments. We will refer to this observation 
in Sect. [6] 

Our 15M runs with varied resolution confirm the trends 
seen for the 11.2 M progenitor. For a neutrino luminosity 



of L v 



2.1 



10 52 ergs _1 , 



for which neither 2D nor 3D sim- 



ulations with standard resolution produce an explosion, we 
find that 2D models with angular binning of 1.5° or better do 
explode, whereas explosions in 3D cannot be obtained with 
angular zones in the range from 1.5° to 3° (Fig. [7J middle 
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Fig. 10. — Upper row: Quasi-three-dimensional visualization of the 11.2M Q simulations in 2D (upper left panel) and 3D (upper right panel) with an electron- 
neutrino luminosity of L Ve = 1.0 ■ 10 52 ergs -1 and an angular resolution of 2°, comparing the structure at 700ms p.b., roughly 150ms after the onset of the 
explosions. Since the explosion started slightly earlier in the 2D model (see the upper panel of Fig.|7]and Table|2j the shock is more extended in the left image. 
While in this case the shock possesses a much stronger dipolar deformation component than in 3D (cf. Fig. [9] lower right panel), the distribution of accretion 
funnels and plumes of neutrino-heated matter exhibits a hemispheric asymmetry in both cases. Because of the axisymmetry of the 2D geometry this concerns the 
hemispheres above and below the ,v-y-plane in the upper left plot, whereas the virtual equator lies in the plane connecting the upper left and lower right corners of 
the top right image and the lower left and upper right corners of the bottom right picture. Note that the jet-like axis feature in the upper left figure is a consequence 
of the symmetry constraints of the 2D setup, which redirect flows moving towards the polar grid axis. Such artifacts do not occur in the 3D simulation despite the 
use of a polar coordinate grid there, too. Lower row: Ray-tracing and volume-rendering images of the three-dimensional explosion of the 11.2 M Q progenitor for 
the same simulation and time displayed in the upper right image. The left lower panel visualizes the outer boundaries of the buoyant bubbles of neutrino-heated 
gas and the outward driven shock, which can be recognized as a nearly transparent, enveloping surface. The visualization uses the fact that both are entropy 
discontinuities in the flow. The infalling matter in the preshock region appears as diffuse, nebular cloud. The right lower panel displays the interior structure by 
the entropy per nucleon of the plasma (red, yellow, green, light blue, dark blue correspond to decreasing values) within the volume formed by the high-entropy 
bubbles, whose surface is cut open by removing a wide cone facing the observer. Note the clear dipolar anisotropy with stronger explosion towards the north-west 
direction and more accretion at the south-east side of the structure. 



panel; Table [2]). The 2D simulations exhibit violent SASI 
sloshing motions and the quasi-periodic appearance of large 
shock asymmetries (Fig. [8] middle panel), and the 2D model 
with 1.5° angular zoning explodes with a huge prolate defor- 
mation (Fig. ITTb . A similar behavior is seen for the 15 M 



runs with L v 



2.2 



10 52 ergs While all 2D models com- 



puted with angular zone sizes between 0.5° and 3° explode, 
we observe an explosion for the 3D calculation with 3° but 
none for the case with 2° angular binning (Tableland Fig. [21 



bottom panel). It is highly interesting that the 3° case, where 
the explosion occurs more readily in 3D than in 2D, is as- 
sociated with a large asphericity of the supernova shock at 
the time the 3D run begins to develop the successful blast 
(Fig. [8j bottom panel). Note again that the structures of the 
higher-resolved 3D model in Fig. [TT] reveal finer details and 
fragmentation into smaller filaments than the corresponding 
2D simulation, despite both having the same zone sizes in the 
angular directions. 
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Fig. 11. — Snapshots of the post-bounce evolution of the 15M G model with an electron-neutrino luminosity of L Ve = 2.1 ■ 10 52 ergs~' and angular resolution 
of 1.5° at fpb = 250, 350, 550, and 700 ms. The color coding represents the entropy per nucleon of the stellar gas. The left half of each panel shows the entropy 
distribution for a 2D simulation, and the right half displays a meridional cut from the corresponding 3D model. While the 2D ran with the given resolution leads 
to a highly prolate explosion, the 3D calculation does not end in a successful blast (see Fig. [7] middle panel, and Table[2). Note that the convective plumes are 
considerably smaller and more fine-structured in the 3D simulation. 



The data listed in Table [2] contain the clear message that 
2D models with better angular resolution usually develop ex- 
plosions earlier in contrast to 3D runs, which explode later 
or not at all when the angular zoning is finer. There can be 
2D exceptions to the general trend (e.g., the 15 M Q cases with 
L Ve = 2.1 ■ 10 52 erg s" 1 and 0.5° and 1° resolution for 400 ra- 
dial zones; see also Fig. [7] middle panel), which are either 
affected by the difficulty to exactly determine the onset of the 
blast in cases with a highly deformed shock, or which can be 
stochastic outliers associated with the chaotic processes lead- 
ing to the explosion. Note in this context that at late times 
M(t) is very flat and therefore differences in f exp correspond 
to only small differences in M exp . For this reason the begin- 
ning of the explosion can be shifted by minor perturbations, 
e.g. connected to stochastic fluctuations. It is also possible 
that for special circumstances the symmetry axis of the 2D 
geometry has an influence on such a non-standard behavior 
because of its effect to redirect converging flows outwards or 
inwards and thus to have a positive feedback on the violence 
of the SASI activity. 

Improved radial resolution for fixed angular grid turns out 
to have a negative influence on the possibility of an explosion 
also in multi-dimensional simulations. 2D runs with better 
radial zoning (600 or even 800 instead of 400 radial zones) 



fail to develop explosions or explode significantly later than 
their less well resolved counterparts (see the 15 M results for 
L Ve = 2.1 -10 52 ergs" 1 and 2.2 -10 52 erg s" 1 in Table©. In gen- 
eral, in ID, 2D, and 3D improved radial resolution shifts the 
onset of the explosion to later times monotonically. The only 
successful 3D model in the set that is useful for the present 
discussion, a 15 M Q run with L Vr = 2.2 ■ 10 52 ergs~' and 3° 
angular resolution, supports our findings in 2D. It shows an 
increasingly delayed explosion for better radial resolution, al- 
though the results with 600 and 800 radial zones appear to be 
nearly converged. All non-exploding 3D models do not yield 
successes also with higher radial resolution. 

The conclusion that good radial resolution is very impor- 
tant for reliable results, in particular when the explosion 
is "ma rginal", would not be a surprise, because ISato et all 
(120091) have pointed out the importance of the radial zoning 
close to the neutron star and around the supernova shock in or- 
der to accurately capture the entropy and vorticity production 
at the shock and to determine growth times and oscillation fre- 
quencies of the SASI. The latter is unquestionably an essential 
ingredient for the success of the neutrino-driven mechanism 
in our 2D runs and it may as well be a crucial component 
for the mechanism to work in 3D. As mentioned in Sect. [H 
however, the sensitive influence of the radial zoning in the 
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discussed model set is mainly a consequence of an artificial 
density peak developing in the neutrino-loss region because of 
the use of the simplified neutrino-cooling treatment. The nar- 
row shape of this numerical artifact in the density structure, 
which enhances the energy emission by neutrinos, can even 
cause a dependence of the results on the particular choice of 
the grid-cell locations. Different from the ID runs the results 
of multi-dimensional simulations with more than ~600 radial 
zones do not seem to converge. Since the local density max- 
imum lies between 10 12 gcirT 3 and 10 13 gcirT 3 in the core 
that is treated spherically symmetrically in our simulations, 
we interpret this phenomenon as the consequence of a sub- 
tle feedback between higher zoning and cooling strength on 
the one hand and multi-dimensional processes in the accre- 
tion layer on the other hand. Because of the artificial nature 
of the underlying density feature, however, we did not further 
explore this finding. 

Fina lly, we remark that prior to our present work IScheckl 
d2007l) has already performed resolution studies with a large 
set of 2D simulations, in which he varied the lateral zone 
width between 0.5° and 4°. In addition, he conducted three 
3D simulations with angular bin sizes of 2—4°. However, in- 
stead of the highly simplified heating and cooling descrip- 
tion used by us he employed the much more sophisticated 
ap proximation for grey neutrino transport described in detail 
in IScheck et al.l ( 120061) . This approximation included, e.g., 
the feedback of accretion on the neutrino emission properties 
and on the corresponding energy and lepton number trans- 
port by neutrinos and antineutrinos of all flavors, as well as a 
more elaborate description of neutrino-matter interactions in 
detailed de pendence on the thermodynamical state of the stel- 
lar plasma. IScheckl (120071) was not interested in a systematic 
exploration of the critical explosion condition, but his project 
was focussed on investigating the possibility of hydrodynamic 
pulsar kicks by successful asymmetric explosions. Despite 
the grave differences of the neutrin o treatments a nd numer- 
ical setups, the results obtained by Scheck ( 2007]) are com- 
patible with our present findings: 2D simulations with higher 
resolution turned out to yield explosions significantly earlier 
and thus also more energetically than the low- resolution runs . 
Within the tested range of angular resolutions IScheckl (120071) 
did not observe any significant differences between his 3D 
models. This, however, may just be a consequence of the fact 
that the models were clearly above the threshold conditions 
for an explosion and did not linger along the borderline be- 
tween blast and failure. 

6. INTERPRETATION AND DISCUSSION 

In this section we discuss the meaning of our results in com- 
parison to previous studies and present an interpretation that 
could explain the main trends found in our multi-dimensional 
simulations with varied resolution. 

6.1. Variation with dimension 

In Sects. |4]and|5] we have repo rted that our simulation s do 
not support the central finding bv lNordhaus et al.l (120101) that 
the tendency to explode is a monotonically increasing func- 
tion of dimension. While we confirm more favorable explo- 
sion conditions in 2D than in ID, we do not observe that in 3D 
considerably lower driving luminosities are needed for a suc- 
cess of the neutrino-driven mech anism than in 2D . More over, 
we cannot confirm the finding bv lNordhaus et al.l d2010l) that 
the mass-weighted average of the entropy per nucleon in the 
gain region, (s(t)), is a quantity that is suitable as an indicator 
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Fig. 12. — Time evolution of the mass in the gain region (in seconds after 
bounce) for simulations in ID (thin dotted line), 2D (thin solid lines), and 3D 
(thick lines). The multi-dimensional models are displayed for all employed 
angular resolutions depicted by different colors. The top panel shows the 
results for the 11.2 M B star with an electron-neutrino luminosity of L,, f = 

1.0 ■ 10 52 erg s _I , the middle panel the results for the 15 M B runs with L Ve = 

2.1 ■ 10 52 < 
10 5: 



■ergs , and the bottom panel the 15A/ Q models for L Ve = 2.2 



n52 er g s 1 -p ne djff eren t cases are the same as in Figs.|7jand[8] 



of the proximity of models to an explosion and thus can serve 
as an explanation of differences between ID, 2D, and 3D sim- 
ulations. In particular, our 3D models turned out to have 
slightly higher mean entropies than corresponding 2D cases 
(Fig.|?li without developing better explosion conditions. This 
raises the question why our models, and multi-dimensional 
simulations in general, have produced successful explosions 
by the neutrino-heating mechanism when corresponding ID 
models fail? 

It is by no means obvious that (s(t)) should increase in the 
gain layer in the multi-dimensional case. While neutrino en- 
ergy deposition naturally leads to a rise of the entropy of the 
heated gas, the averaging process over the volume of the gain 
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Fig. 13. — Analogous to Fig. [12] but for the time evolution of the total net 
rate of neutrino heating in the gain region. 

layer also encompasses the downdrafts carrying cool matter 
from the postshock region towards the gain radius and the 
neutron star. These downdrafts are much denser, they are 
hardly heated by neutrinos because of their extremely rapid 
infall, and they can contain more mass than the surrounding, 
dilute bubbles that are inflated by the expanding, neutrino- 
heated plasma. It is therefore not clear that the spatial (mass- 
weighted) average (s{t)) grows in multi-dimensions compared 
to ID runs. 

Moreover, it is not even clear that convective overturn in the 
gain layer must lead to an average entropy of the neutrino- 
heated gas itself that is higher than in ID simulations. Dif- 
ferent from the ID case high-entropy matter becomes buoy- 
ant and begins to float in the multi-dimensional environment. 
Thus the heated gas is quickly carried away from the vicinity 
of the gain radius, where neutrino-energy deposition is maxi- 
mal, to larger radii. Such dynamics of the gas can well limit 
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Fig. 14. — Analogous to Fig. 1121 but for the time evolution of the ratio of 
advection timescale to heating timescale in the gain layer. 

the amount of energy and entropy that is stored in individual 
chunks of matter. Little, if any of the gas is subject to mul- 
tiple overturn cycles bringing the gas close to the gain radius 
more t han once as suggeste d by the "convecti ve engine" pic- 
ture of iHerant et aT] (119941) but questioned by iBurrows et al.l 
( 1 19951) . Instead, the majority of the heated gas either expands 
to larger distances, pushing shock expansion, or, in the disad- 
vantageous case, is swept back below the gain radius (e.g. by 
large-amplitude sloshing motions of the shock), where it loses 
its energy again by efficiently reradiating neutrinofl 
Our results imply that the dominant effect that makes the 

3 The real multi-dimensional situation is even more complicated. The 
mentioned violent sloshing motions of the shock can cause strong shock - 
heating of the postsh ock m atter as discussed in detai l by[S check et alj 120081) , 
IBlondin et al. (2003), and Blondin & Mezzacapplj )2006bf) . thus not only 
massively affecting the postshock flow but also providing an additional en- 
tropy source besides neutrino-energy deposition. 
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Fig. 15. — Kinetic energies of angular mass motions in the gain layer as 
functions of time after bounce for the 11.2 M Q runs with an electron-neutrino 
luminosity of L Ve = 1.0 ■ 10 52 ergs -1 (top panel) and the 15 M Q runs with 
L Ve = 2.1 ■ 10 52 ergs~' (middle panel) and L„ f = 2.2- lO^ergs -1 (bottom 
panel). Thin solid lines correspond to the lateral kinetic energy of 2D mod- 
els, while for 3D simulations (thick lines) the lateral, azimuthal, and total 
kinetic energies are represented by dotted, dashed, and solid line styles, re- 
spectively. Both angular directions contribute essentially equally to the t otal 
kin etic energy of nonradial motions in the 3D case. As in Figs. 171181 [HI 
and !13l different colors depict different angular resolutions. It is visible that 
for models closer to a success of the neutrino-driven mechanism the angu- 
lar kinetic energy exhibits larger temporal variations and an overall trend of 
increase as the onset of the explosion is approached. 

multi-dimensional case more favorable for an explosion than 
spherical symmetry is associated with an inflation of the 
shock radius and postshock layer, driven by the buoyant rise 
and expansion of the plumes of neutrino-heated plasma. In 
course of the postshock volume becoming more extended, the 
integrated mass M ga j n in the heating layer increases compared 
to the ID case. This can be seen in Fig. [121 which displays 
the mass in the gain layer as function of post-bounce time for 



the 1 1.2 and 15 M runs with the different neutrino luminosi- 
ties and resolutions already shown in Figs.[7]and[8] While the 
mass-averaged entropy (s) in the gain region hardly changes, 
the integral value of the entropy, M ga j n (.y), clearly increases 
with models coming closer to explosion. This dependence is 
particularly well visible when 2D and 3D models with differ- 
ent resolutions are compared with each other. 

The longer dwell times of matter in the gain l ayer o f 
2D simulations observed by iMurphv & Burrows! d2008l) . 
which correspond t o th e advection times T a d v evaluated by 
iBuras et alj d2006bl) and iMarek & Jankal j2009l) (though dif- 
ferent ways of calculation have been considered, in particu- 
lar for the multi-dimensional case), are a manifestation that 
a growing mass accumulates in the postshock region to get 
energy-loaded by neutrino absorption and to finally drive the 
successful supernova blast. In near-steady-state conditions 
the mass accretion rate through the gain layer is equal to 
the mass infall rate M ahead of the shock, where it is de- 
termined by the cor e structure of the prog enitor star. Since 
T adv ~ Mg a i n /M (cf. IMarek & Jankal 120091) a larger value of 
r a dv correlates with a higher mass M aa ; n . Accordingly, the 
total net heating rate Qgain and thus the heating efficiency 
e = Q g ain/(L Ve + Ly e ) = <2 g ain/(2L v ) of the gas residing in the 
gain layer is als o higher for models that de velop an explosion 
(see Fig. [131 and IMurphv & Burrowsll2008l) . Figure [14] shows 
the ratio of the advection (dwell) timescale to the neutrino- 
heating timescale in the gain layer (our evalua tion employs 
the N ewtonian analog of the formulas given in iMiiller et alJ 
120121) for the models also displayed in Figs. fl2l and [131 The 
same trends as in the previous images can be seen. Models 
closer to an explosion exhibit higher values of the timescale 
ratio. The ratio approaches unity, indicating the proximity to a 
runaway instability, roughly around the time when we define 
the onset of an explosion (i.e., when the average shock ra- 
dius of Fig.[7]passes 400 km). The quality of the coincidence 
of these moments, however, differs from model to model and 
depends on the exact definition used for the timescales and the 
accuracy at which the relevant quantities can be evaluated in 
highly perturbed flows. Nevertheless, the evolutionary trend 
of the timescale ratio (increasing or decreasing) is suggestive 
for whether a simulation run leads to a successful explosion 
or failure. 

A larger mass in the gain layer and higher total net energy 
deposition rate are therefore better indicators of the proxim- 
ity of our models to explosion than the mean entropy of the 
gas in this region, which does not exhibit th e 1D-2D-3D hi- 
erarch y with dimension found previously by iNordhaus et al. I 
(120101) . As discussed in Sect. [4] the main reason for this 
discrepancy are most probably the different treatments of 
neutrino lepton number losses and our consequential recal- 
ibration of the energy source terms. This leads to signifi- 
cantly higher energy drain from the cooling layer in our sim- 
ulations. While this hypothesis is supported by tests that 
we conducted in ID, we cannot be absolutely certain that 
no other effects play a role for the discrepan cies between 
our results and those of INordhaus et al.l (120101) . because de- 
tailed cross-comparisons are not available and our knowledge 
of the details o f the i mplementation of neutrino effects by 
INordhaus et al~l (120101) may be incomplete. Other potential 
reasons for differences may be connected to the hydrodynam- 
ics scheme (Prometheus with a higher-order Godunov solver 
and directional splitti ng vs. Castro with unsplit methodology, 
lAlmgren et al.l I20T0I) . the employed grid (polar coordinates 
vs. structured grid with adaptive refinement by a nested hi- 
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erarchy of rectangular grids), potentially — though not very 
likely — the use of a ID core above 10 12 g cirT 3 in our simula- 
tions, or differences in the exact structure and properties of the 
infall region upstream of the stalled shock as a consequence 
of different treatments of the collapse phase until 15 ms af- 
ter core bounce (due to full neutrino transport plus a nuclear- 
burning approximation in the Prome theus- Vertex code vs. 
the simple deleptonization scheme of Liebendorfer 2 0051) or 
of different seeding of nonradial hydrodynamic instabilities 
(in our case by imposed, small random seed perturbations of 
the density), or linked to differences of th e low-dens it y EoS 
outside of the application regime of the LShen et al.l (U998I) 
EoS. 

Despite these uncertainties about the exact cause of the 
differences, whose ultimate elimination will require system- 
atic and time-consuming studies, our results, as they are, 
send a clear message: The outcome of the 1D-2D-3D com- 
parison and the effects of the third dimension advertised by 
iNordhaus etaTI(l2010h ' 'as a key to the neutrino mechanism 
of core-collapse supernova explosions" are not at all robust 
results. Instead, the exact slope of the critical explosion con- 
dition L V {M), its location, and its shift with dimension, as 
well as the existence of a 1D-2D-3D hierarchy of the mass- 
averaged entropy in the gain layer seem to depend sensitively 
on subtle details of the neutrino treatment or other numerical 
aspects of the simulations. 

6.2. Resolution dependence 

Let us now turn to the second, highly interesting question 
in connection with our set of models, namely to the resolu- 
tion dependence of our results. Our set of simulations per- 
formed with different angular binnings reveals that quantities 
that turned out to diagnose healthy conditions for an explo- 
sion, i.e. the growth of the average shock radius, the degree of 
shock deformation, or the mass and total heating in the gain 
layer (but not the mass-averaged entropy of the matter in the 
gain region), show a clear dependence on the angular zoning 
(see Figs. [7] [8] [12] and[13]in contrast to Fig.|4|. In particular, 
2D models with better angular resolution exhibit more favor- 
able conditions and explode more readily (in agreement with 
results obtained by IScheckl 120071 with a more sophisticated 
treatment of neutrino transport than the simple heating and 
cooling source terms applied in our investigation), whereas 
3D models obey the opposite behavior. What does the reso- 
lution dependence of our simulations tell us about the mech- 
anism leading to explosions in our models? And how can we 
understand the puzzling finding that 2D and 3D runs follow 
opposite trends when the angular resolution is refined? 

We interpret this as a manifestation of two aspects or facts: 

(1) The success of our models, at least in the neighborhood 
of the explosion threshold, is fostered mainly by large- 
scale mass flows as associated with strong SASI activ- 
ity, but not by enhanced fragmentation of structures and 
vortex motions on small spatial scales. 

(2) Our resolution study reflects the consequences of the 
turbulent energy cascade, which redistributes energy 
fed into the flow by external sources in opposite direc- 
tions in 2D and 3D: While in 3D the turbulent energy 
flow goes from large to small scales, it pumps energy 
from small to large spatial scales in 2D. 

Point (1) is supported by the kinetic energies of nonradial 
mass motions in the gain layer of the 2D and 3D models plot- 



ted in Fig. [15] From this picture it is obvious that in the case 
of successful models the angular kinetic energy is higher and 
shows an overall trend of growth in time until the blast has 
taken off. Moreover, the spiky maxima and minima of quasi- 
periodic variations, which are indicative of the presence of 
low-order SASI modes, are significantly larger for exploding 
models. This does not only hold for 2D models, whose lateral 
kinetic energies exhibit variation amplitudes of several 10% 
and partially up to even ~50% of the time-averaged value. It 
is also true for 3D models, although in this case the ampli- 
tudes are generally smaller and the nonradial kinetic energy 
is split essentially equally into lateral and azimuthal contribu- 
tions. When comparing successful runs in 2D with those in 
3D, our studies suggest that in both cases the shock exhibits a 
growing degree of asphericity (expressed by the standard de- 
viation of the shock deformation plotted in Fig. [8]) when the 
explosion is approached, and the kinetic energy of nonradial 
mass motions reaches roughly the same magnitude (Fig. [15]), 
at least for models near the explosion threshold. 

Actually a variety of observations can be interpreted as sup- 
port of the hypothesis that flows on the largest possible scales 
rather than on small scales play a crucial role for the success 
of the neutrino-heating mechanism in our simulations: 

• The strength of low-mode SASI activity in 2D models 
as indicated by growing fluctuations of the angular ki- 
netic energy and of the shock deformation (Figs. Q3] 
[8]) increases with higher resolution in clear correlation 
with an earlier onset of the explosion (Fig. [7). Stronger 
SASI activity obviously facilitates explosions, which is 
visible by a growing average shock radius as well as 
larger mass and higher total heating rate in the gain 
layer. 

• Exploding models in 2D as well as in 3D exhibit large 
shock deformation at the time of explosion (although 
the relative asphericity os/Rs of the shock surface is 
somewhat smaller in 3D than in 2D; see Figs.[7]and[8]). 

• More fine structure on small spatial scales, which can 
be seen in 3D models computed with higher resolution 
in Figs. |9hTT1 does not imply improved conditions for 
an explosion. 

• Exploding 2D models are not connected with the high- 
est mean entropies in the gain region (Fig. [4]). This fuels 
doubts in a random-walk picture where turbulent vortex 
motions on small scale s enlarge the residence tim e of 
matter in the gain layer dMurphv & Burrowsll2008l) and 
thus could allow for more energy absorption of such 
mass elements from neutrinos. If this effect occurred, 
it does not concern a major fraction of the mass in the 
gain region. 

Point (2) is the only plausible argument we can give for 
explaining the opposite response to higher angular resolution 
that we discovered in our 2D and 3D simulations. The se- 
quence of 2D runs with gradually reduced lateral zone sizes 
reflects the growing violence of large-scale flows by higher 
fluctuation amplitudes of the kinetic energy in the gain layer 
(Fig.fTBIl and larger temporal variations of the average shock 
radius and shock deformation (Figs. [7] [8). In contrast, more 
energy on small spatial scales in the 3D case manifests itself 
by a progressing fragmentation of the flow, leading to a grow- 
ing richness of vortex structures and finer filaments in the case 
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Fig. 16. — Turbulent energy spectra E(l) as functions of the multipole order / for different angular resolution. The spectra are based on a decomposition of 
the azimuthal velocity vg into spherical harmonics at radius r = 150 km and 400 ms post-bounce time for 15M Q runs with an electron-neutrino luminosity of 
L Ve = 2.2 ■ 10 52 ergs~'. Left: 2D models with different angular resolution (black, different thickness) and, for comparison, the 3D model with the highest 
employed angular resolution (grey). Right: 3D models with different angular resolution and, for comparison, the 2D model with the highest employed angular 
resolution (grey). The power-law dependence and direction of the energy and enstrophy cascades (see text) are indicated by red lines and labels for 2D in the 
left panel and 3D in the right panel. The left vertical, dotted line roughly marks the energy-injection scale, and the right vertical, dotted line denotes the onset of 
dissipation at high / for the best displayed resolution. 

indeed confirm the predictions from 3D and (planar) 2D tur- 
bulence theory, at least for sufficiently high multipole order 



of 3D models with smaller angle bins (Figs.l9l-[TT1>. As a con- 
sequence, 3D models with higher angular resolution become 
more similar to the ID case in various quantities that we con- 
sidered as explosion indicators, see, e.g., Figs. [7j [8j [T2l and 

Both the powerful coherent mass motions of the S ASI layer 
in 2D and the vivid activity in small vortex structures in the 
3D environment are fed by two external sources which sup- 
ply the postshock layer with an inflow of fresh energy: (i) 
gravitational potential energy that is released by the continu- 
ous stream of matter falling through the accretion shock and 
(ii) energy deposition by neutrinos. The energy stored in the 
fluid is then redistributed towards small or large scales ac- 
cording to the turbulent cascades characteristic of two- and 
three-dimensional environments. 

Direct evidence for the action of different turbulent energy 
cascades in 2D and 3D can be obtained by considering the 
energy spectrum E(k) of turbulent motions as a function of 
wavenumber k in the gain region. The spectral shape of E(k) 
can already be adequately established by considering only the 
azimuthal velocity vg at a given radius using a decomposition 
into spherical harmonics Yi m (0,<p): 



i 

I 

m=-l 



f Kn 

Jn 



(6,<f>) yfpvg(r,6,4>)dQ, 



(8) 



Here, the velocity fluctuations have been expressed in terms 
of the multipole order / instead of the wave number k. A sum- 
mation over the energies of modes with the same I has been 
carried oufl, and in order to obtain smoother spectra, we av- 
erage E{1) over 30 km in radius and over 10 timesteps. One 
expects that the resulting spectrum E(l) directly reflects the 
properties of E(k) such as the slopes in different regimes of 
the turbulent cascade^. The computed spectra E(l) (Fig. [Tol l 

4 Note that a factor y[p has been introduced to ensure that the integrated 
energy of all modes sums up to the total kinetic energy contained in az- 
imuthal motions at radius r = 150 km (modulo a normalization factor) (cp. 
lEndeve et aU2012t) . 

5 For the precise relation between Fourier and spherical ha rmonics power 
spectra, the reader may consult Chapter 21 of Peebles ( 1993). For a power- 
law spectrum E(k) oc k a , one obtains E(l) cc (21 + l)T(l + a/2 + l/2)/T(/ - 
a/2 + 3/2), or E(l) oc l" in the limit of large /. In practice, the power-law 
indices of E(I) and E(k) appear to correspond well to each other already for 



I. 

In 3D (right panel), a p ower-law spectrum with E(l) cc l~ 5 ' 3 
dLandau & Lifshitzl 19591) develops at intermediate wavenum- 
bers as the resolution is increased, reflecting the transfer of 
energy from large to small scales in a forward cascade un- 
til dissipation takes over at large /. At high resolution, the 
energy contained in small-scale disturbances increases, as the 
dissipation range moves to larger /. One observes that the 5/3- 
power-law is broken at low / (I < 10), suggesting that kinetic 
energy is injected into the flow at wavenumbers I as 10, i.e. at 
scales typical for growing convective plumes. 

By contrast, the power-law dependence E(l) oc r 5/3 approx- 
imately holds for I < 10 in 2D as a result of the reverse en- 
ergy cascade (iKraichnanl 1 967b . The energy injected at I * 10 
is therefore not transferred to the dissipative range; only en- 
strophy (the squared vorticity of the velocity field) is trans- 
ported in a forward cascade with a different power-law index 

(£(0 oc r 3 ). 

This appears to be a natural explanation for the predomi- 
nance of large-scale and small-scale structures in 2D and 3D, 
respectively. Moreover, this picture suggests that as dissipa- 
tion affects the "injection scale" at / as 10 less with increasing 
resolution, the differences between 2D and 3D become more 
pronounced with finer grid zoning. 

Our analysis of the spectral properties of turbulence thus 
further strengthens our view that the trends seen in our sim- 
ulations strongly suggest that nonradial kinetic energy avail- 
able on large scales, not on small scales, assists the devel- 
opment of an explosion by the neutrino-heating mechanism. 
This explains why 2D models with higher angular resolution 
tend to explode earlier and thus at higher values of the mass- 
accretion rate than less resolved models. On the other hand, 
the energy "drain" by vortex motions on ever smaller scales 
— with the same reservoir of pumping energy per unit mass 
being available from accretion and neutrino heating — disfa- 
vors explosions in better resolved 3D models. 

We therefore conclude that the key to the mechanism of 
core-collapse supernova explosions seems intrinsically and 
tightly linked to the question how much kinetic energy of 

/ > 4. Empirically, broken power laws transform in a similar manner. 
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the matter in the gain region can be accumulated in nonra- 
dial fluid motions on the largest possible scales, i.e., in the 
lowest-order spherical harmonics modes of nonradial hydro- 
dynamic instabilities. The predominant growth of such flows 
is typical of SASI activity, whose lowest-order spherical har- 
monic s modes possess the highest growth rates dBlondin et al J 
20031; iBlondin & Mezzacappal |2006bi iFoglizzo et al.l 120061 



20071 lOhnishi et alj|2006h . Strong SASI motions drive shock 
expansion, increase the gain layer and its mass content, al- 
low a larger fraction of the accreted matter to stay in the 
gain layer and be exposed to efficient ne utrino heating, and 
thus aid the developm ent of an explosion (IScheck et al.ll2008t 
iMarek & Jankal 120091) . However, our models do not show a 
systematic trend of higher average entropies of the matter in 
the gain layer for models closer to explosion. Instead, we find 
that such models have larger mass, larger nonradial kinetic 
energy, larger total neutrino-heating rate, and larger total en- 
tropy in the gain layer. 

7. SUMMARY AND CONCLUSIONS 

We have performed a systematic study of the post-bounce 
evolution of supernova cores of 1 1.2 and 15 M Q and their ex- 
plosion by the neutrino-heating mechanism in ID, 2D, and 
3D, employing simple neutrino cooling and heating terms 
with varied values of the driving luminosity. We conceptu- 
ally followed previous stud ies by lMurphy & Burrows! (120081) 
and lNordhaus et al~1 (120101) . but did not apply the deleptoniza - 
tion treatment that they adopted from iLiebendorfefi d2005l) . 
who introduced it for an approximative description of neu- 
trino losses during the infall phase until core bounce. We ar- 
gued (Sect. |2| that this a pproximation —with or without the 
source term proposed by iLiebendorferl (120051) to account for 
entropy generation in neutrino-electron scatterings — does not 
provide a suitable treatment of the evolution of the electron 
abundance after core bounce. Therefore we did not consider 
changes of the net electron fraction Y e of the stellar plasma 
at times later than 15 ms after bounce, up to which the col- 
lapse was followed with the Prometheus-Vertex code includ- 
ing full neutrino transport. While ignoring Y e changes subse- 
quently is certainly not a good approximation, it is not neces- 
sarily more unrealistic than describing the lepton-number evo- 
lution dur ing the accretion pha se of the stalled shock by the 
scheme of ILiebendorferl J2005). As a consequence, we had to 
replace an expone ntial factor e~ T e S , whi ch w as introduced in 
an ad hoc way by iMurphy et al.1 (120091) and iNordhaus et alTI 
(120101) to damp the neutrino source terms at high optical 
depths r e ff, by e' Tclr ^ 2J in order to reproduce the minimum 
luminosity found to yield e xplo sions in the ID simulati ons by 
IMurphy & Burrows! (120081) and lNordhaus etal.l (120101) . This 
modification led to enhanced neutrino losses in the cooling 
layer, which were better compatible with total energy loss 
rates fo und in simulations with detailed neutrino transport, 
e.g., in iBuras et ail d2006al) . and is responsible for some of 
the findings and differences discussed in Sect. [4] 

Our results and conclusions can be briefly summarized as 
follows: 

1 . We cannot reproduce the exact slopes and relative loca- 
tions of the critical c urves L V (M) of ID, 2D , and 3D 
simulations found by INordhaus et al. I (120101) . While 
our results confirm the well-known fact that explosions 
in 2D occur for a lower driving luminosity L v than in 
ID when the mass accretion rate M is fixed, we cannot 
discover any significant further reduction when we go 



from 2D to 3D. 

2. We cannot confirm that the mass-averaged entropy of 
the matter in the gain region, (s), is a good diagnostic 
quantity for t he p roximity to an explosion. As we ar- 
gued in Sect. 16.11 it is neither clear nor necessary that 
(s) is higher for cases where explosions are obtained 
more readily. Our successful 2D models do not exhibit 
larger mean entropies than the corresponding ID cases, 
which fail to explode. Instead, we observed that the to- 
tal mass, total entropy, total neutrino-heating rate, and 
the nonradial kinetic energy in the gain layer are higher 
in cases that develop an explosion. 

3. We conclude that the tendency for an explosion as 
a monotonically increasing function of dimension as 
well as the 1D - 2D-3D hierarchy of (s(t)) found by 
INordhaus et al. I (120101) are not robust results. They 
seem to be sensitive to subtle differences of the ap- 
proximations of neutrino effects (and/or to other differ- 
ences in the numerical treatments of the models). It 
is therefore unclear how far studies with radical sim- 
plifications of the neutrino physics (without detailed 
energy and lepton-number source terms and transport; 
no feedback between accretion and neutrino proper- 
ties) can yield results that are finally conclusive for the 
explosion-triggering processes in real supernova cores. 

4. Increasing the angular resolution we observed a clear 
tendency of 2D models to e xplode earlier, in agree- 
ment with previous results bv lScheckl (120071) . who em- 
ployed a more sophisticated treatment of neutrino ef- 
fe cts based on t he tra nsport approximation described 
in IScheck et al.l (120061) . In contrast, 3D models show 
the opposite trend and in a variety of quantities and 
aspects become more similar to their ID counterparts. 
The easier explosion of the 2D models is connected to 
an enhanced violence of large-scale mass motions in 
the postshock region due to SASI activity, whereas 3D 
models with better angular resolution appear to develop 
less strength in low-order SASI modes. 

5. We interpret this finding as a consequence of the op- 
posite turbulent energy cascades in 2D and 3D. In 2D 
the energy continuously pumped into the gain layer by 
neutrino heating and the release of gravitational binding 
energy flows from small to large scales and thus helps 
to power coherent mass motions on the largest possible 
spatial scales. In contrast, in 3D this energy seems to 
instigate flow vorticity and fragmentation of structures 
on small scales. Evidence for this interpretation is pro- 
vided by Fig.fTSI 

6. We also conclude from our resolution studies that the 
presence of violent mass motions connected to low- 
order SASI modes is favorable for a n explosion (in 
agree me nt with arguments given by IMarek & Jankal 
120091 and IScheck et alj|2008h . This is supported by the 
fact that 2D and 3D models that are closer to explosion 
show signs of growing power in large-scale mass mo- 
tions (signalled by growing fluctuations of the kinetic 
energy of nonradial velocity components) and in partic- 
ular develop significant shock deformation and global 
ejecta asymmetries when the explosion sets in. 
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7. We found that higher radial resolution makes explo- 
sions more difficult with the setup chosen for the inves- 
tigated set of models. Higher resolution turned out to 
prevent explosions or to let them occur later in simula- 
tions in ID, 2D, and 3D. This result could be diagnosed 
to be a consequence of a local density maximum in the 
neutrino-cooling layer, which grows with higher reso- 
lution and enhances the energy loss by neutrino emis- 
sion. This density peak, however, is a numerical artifact 
of the employed simple neutrino-cooling treatment by 
an analytic source term which is exponentially damped 
at high densities. 

The lack of very precise information on the physics in- 
gredients and their exact implementation, e.g., details of the 
treatment of neutrino source terms, low-density EoS, and pro- 
genitor data when mapped into the simulation and seeded 
with small random perturbations, as well as a variety of 
methodical differences like the hydrodynamics scheme, nu- 
merical grid, and the use of a ID core at high densities or 
not, prevent us from presenting a rigorous proof that could 
ca usally link the d i screpa ncies between our results and those 
o flNord haus et aT~l (120 1 Ot) to one or more well understood rea- 
sons. We think that the nagging uncertainties in this context 
demand a future, involved, collaborative code-comparison 
project. This will also require considerable amounts of com- 
puter time for further 3D simulations, in particular with high 
resolution, thus needing more computer resources than avail- 
able to us for the described project. 

Despite this deficiency, however, our results suggest that 
th e differences of 3D com pared to 2D simulations observed 
bv lNordhaus et al.l (120101) are unlikely to be a robust outcome 
but seem to depend on relevant aspects of the modeling (most 
probably the neutrino physics but potentially, and not finally 
excluded, also technical aspects). 

We therefore conclude that the influence of 3D effects on 
the supernova mechanism is presently not clear. We strongly 
emphasize, however, that the fact that our results do not cor- 
roborate improved explosion conditions in 3D compared to 
2D cannot be used as an argument that 3D effects do not fa- 
cilitate the supernova explosion mechanism or are of minor 
importance. We just think that in the context of the neutrino- 
driven mechanism the relevance and exact role of 3D fluid 
dynamics are not understood y et. We therefore have th e opin- 
ion that the results obtained bv lNordhaus et al.l (120101) do not 
justify their claims that 3D hydrodynamics offers the key to a 
fundamental understanding of the neutrino mechanism while 
other physics in the supernova core, like general relativity 
or the properties of the nuclear EOS, are only of secondary 
importance. Though this may well be right, such statements 
at the present time are premature and not supported by solid 
facts and results. 

Our study, however, raises further important questions. 
How far can our understanding be developed on grounds of 
modeling approaches that employ radical simplifications of 
the neutrino physics? Which aspects of the complex inter- 
play between different components of the problem are linked 
to the essence of how the explosion is triggered by the combi- 
nation of neutrino energy supply and nonradial hydrodynamic 
instabilities? Examples for such mutually dependent compo- 
nents are the neutrino transport and hydrodynamics, the neu- 
tron star core evolution and fluid motions around the neutron 
star, or the mass flux from the accretion shock to the deceler- 
ation layer (both being the coupling regions for the advective- 



acoustic cycl e that is thought to b e responsible for the SASI 
growth; e.g., IScheck et al.l 12008b and the conditions in the 
neutrino heating and cooling layers. Much more work needs 
to be done to find the answers of these questions. 

A major shortcoming of the setup applied in previous works 
and adopted also in our investigation is the neglect of neu- 
trino cooling and deleptonization inside the proto-neutron 
star. The employed simple neutrino source terms and their 
exponential suppression at high optical depths do not al- 
low the neutron star to evolve. Underestimating the neutron 
star contraction, however, slows down the infall velocities in 
the postshock layer and thus has disfavorable co nsequences 
for the growth of the SASI ( IScheck et al J 120081) . similar to 
the effects of reduc ed neutrino losses in the cooling layer 
( IScheck et al.l 120081) or i ncreased nuclear photodissoci ation 
behind the stalled shock dFernandez & Thompsonll2~009l) . On 
the contrary, a slower postshock flow improves the conditions 
for the growth of convective instability, whose development 
is supported by a high ratio (larger than ~3 signals a linearly 
unstable situation) of the advection timescale through the gain 
region divided b y the inverse o f the B r unt- Vaisala f r equenc y 
as discussed by iFoglizzo et all ( 120061) . iBuras et al.l d2006bl) . 
and lScheck et alj ( 120081) . It is therefore well possible that the 
weak S ASI activity diagnose d in recent 2D and 3D simula- 
tions by iBurrows et al.l ( 120121) — leading them to the conclu- 
sion that dipolar asymmetries are caused by convection rather 
than the SASI — is an artifact of the approximative treatment 
of neutrino cooling and the disregard of neutron star contrac- 
tion in their models. In full-sc ale supernova model s with so- 
phisticated neutrino transport, IBuras et ail d2006bl) observed 
differences in the growth conditions for convection compared 
to the SASI in collapsing stellar cores of a variety of progen- 
itor stars. Also the high-density equation of state and general 
relativity, influencing the contraction behavior of the nascent 
neutron star, can make a difference dMarek & Jankal l2009t 
iMiiller et al. 2012). Claims, based on highly simplified mod- 
els, that SASI is less important than convecti on and at most 
a min or feature of the supernova dynamics ( IBurrows et alj 
12011) are therefore certainly premature. 

Finally, our resolution study suggests that the action of the 
turbulent cascade in 3D extracts energy from coherent large- 
scale modes of fluid motion and instead fuels fragmentation 
and enhanced vortex flows on small spatial scales. At least 
in our 3D models with better grid zoning the appearance of 
finer structures in the postshock flow was connected with a 
tendency of damping the development of explosions. While 
a finally convincing proof of such a negative feedback may 
require much better resolved simulations than we presently 
can afford to conduct (in order to minimize numerical dissipa- 
tion on small scales), this result implies that good resolution 
—con siderably higher than recently used by Tak iwaki et alj 
(1201 ll) . whose 3D simulation had only 32 azimuthal zones 
(corresponding to a cell size of 1 1.25°) — is indispensable to 
clarify the 3D effects on the explosion mechanism. More- 
over, our result points to an interesting direction. It is pos- 
sible that the success of the neutrino-driven mechanism in 
3D is tightly coupled to the presence of violent SASI activ- 
ity, a connection that was found before — and is confirmed by 
our present study — to foster explosions in 2D? If so, what 
is the key to instigate such violent SASI motions of the su- 
pernova core in three dimensions? Will they occur with a 
better (more realistic) treatment of the neutrino transport and 
correspondingly altered conditions in the heating and cool- 
ing layers and in the contracting core of the proto-neutron 
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Fig. 17. — Mass accretion rate for the 15 M Q progenitor. The red curve shows the line of Fig. \T\ The black symbols represent the values extracted from our 



simulations at the time i cxp when the explosion sets in. Green symbols are data from Murphy & Burrows (2008) and blue symbols those from Nordhaus et al. 



<201CD . Different symbols are used for results of ID, 2D, and 3D simulations. 

star? Or are they associated with stellar rotation, which even 
with a slow rate can initiate t he faster growth of spiral (non- 
axisymmetric) SASI modes dBlondin & Mezzacappall2006at 
Yamasaki & Foglizzol 120081 : llwakami et al.ll2009b iFernandezI 
20101) ? Or is strong SASI activity in the supernova core trig- 
gered by large- scale inhomogeneities in the three-dimensional 
progenitor star dArnett &~M eakin 201 1), which could provide 
a more efficient seed for SASI growth than the random cell- 
to-cell small-amplitude perturbations employed in our simula- 
tions? Should the presence of large-amplitude SASI mass mo- 
tions indeed turn out to be the key to the neutrino mechanism 
in 3D, it would mean that neutrino-driven explosions are not 
only a generically multi-dimensional phenomenon, but one 
that is generically associated with dominant low-order modes 
of asymmetry and deformation from the very beginning. 

While this paper raises many more questions than it is able 
to answer, it definitely makes clear that our understanding of 
the supernova physics in the third dimension is still in its very 
infancy. A virgin territory with distant horizons lies ahead of 
us and awaits to be explored. 
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APPENDIX 

A. SIMULATIONS WITH PARAMETRIZED DELEPTONIZATION TREATMENT FOR THE CORE-COLLAPSE PHASE 

In th is Appendix we briefly report on our efforts to reproduce the ID results of lMurphv & Burrows! (12008b and lNord haus et al. 
( 120101) for the critical explosion conditions of the 15 M Q progenitor, applying a neutrino treatment that was intended to copy the 
procedure outlined in these publications as closely as possible. 

For this purpose we retained the exponential suppression factor e~ Tdl of Eqs. (0|i and (0 without a reduction factor of 2.7 in 
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Fig. 18. — Trajectories of the electron fraction with density deduced from different core-collapse studies (see text for details) and employed in our ID simulations 
for parametrizing lepton losses by neutrino emission in the stellar core according to Liebendorfer ( 20051) . 

the exponent, and the lepton evolution before and after core bounce was described by employing a predefined Y e (p) relation. We 
also aimed at reproducing the core infall of the previous works as closely as possible, because the density structure of the infall 
region ahead of the stalled shock determines t he mass-infall rate M(t) at the s hock, and some differen ces became visible when 
we compared our values with those given by iMurphv & Burrowsl (120081) and iNordhaus etail (120101) (Fig. Q7}. We therefore 
recomputed the collapse phase from the onset of gravitational instability of the progenitor core through core bounce with th e 
deleptonization scheme of ILiebendorferl (120051) . Entropy changes were taken into account as suggested b y ILiebendorferl (120051) . 
but were switched off after core bounce following lMurphv & Burrowsl (120081 ) and INordhaus et aLl (120101) . 

We tested three different cases for the functional relation Y e {p): First, we used a tabul ated result for the Y e (p) evolu- 
tion as obtained with the Prometheus- Vertex code and state-of-the-art electron-capture rates (lLanganke et al.l 120031) ("Delep- 
tonization 1"). Second, we applied a Y e (p)-table provided by Christian Ott as a co-developer of the CoCoNuT code 
(http://www.mpa-garching.mpg.de/hydro/COCONUT/). These data are based on collaps e simulations with th e Vulcan/2D code 
dLivne et al.ll2.004h ("Deleptonization 2"). Third, we employed a fitting formula given by ILiebendorferl J2005I) for the parameters 
of model G15 ("Deleptonization 3")- All three Y e (p) trajectories are depicted in Fig. [18] 

The three sets of ID simulations conducted for this Appendix were performe d with 800 radial zones. The cor responding criti- 
cal lum inosity curves L Ve (M) are displayed in Fig.[19]in comparison to those of IMurphv & Burrowsl (120081) and INordhaus et af~l 
(120101) and to our results of Fig. [2] (for 400 zone s, because for this r e soluti on the calibration of the exponential suppression 
factor for best agreement with the critical curve o f IMurphv & Burrowsl (120081) was done). The overall slopes of all three cu rves 
are similar but none o f them is quantitatively or qualitatively in good agreement with those of IMurphv & Burrowsl (120081) and 
INordhaus et al. 1(120101) . Explosions in our simulations occurred significantly more readily (i.e., for lower L Ve ) than in the previous 
works. This suggests less cooling in our runs, although we made all possible efforts to exactly follow the description of the neu- 
trino treatment in those papers. The steep rise and partly backward bending of our curves for M values around 0.2-0.3 M Q s _1 can 
be understood by an inspection of Fig. [20] which shows the time evolution of the shock radius for simulations with prescription 
"Deleptonization 1" for a selection of neutrino luminosity values. One can see that in the case of L Ve = 2.95 ■ 10 52 ergs _1 the 
shock makes a larger excursion before it returns again. Its reexpansion, leading to an explosion, therefore happens later than in 
the model with L Ve = 2.7 ■ 10 52 ergs~', where the first shock expansion is much less strong. Correspondingly, the explosion in 
the former case sets in at a later time and lower mass accretion rate than in the latter case, explaining the backward bending of 
L Ve (M) in this regime of luminosities and M. The nearly horizontal parts of the critical curves can be understood by the fact that 
for such high values of the luminosities the neutrino cooling (with the unmodified e~ TcS suppression factor) is so weak that the 
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Fig. 19. — Critical curves for the electron-neutrino luminosity (L Ve ) versus mass accretion rate (M), representing the explosion threshold for different sets 
of ID simulations of the 15M R pro genitor. The black line corresponds to our results shown as black curve in Fig. [2] (see also Table [TJ, red are results of 
IMurphv & Burrows 12008), green of Nordhaus et al. (2010), and the three add i tional curves (dark blue, light blue, and pink) correspond to different sets of 
simula tion s that we performed with the deleptonizati on treatment of Liebendorfer ( 2005) for the core-collapse phase and the different electron-fraction trajectories 
of Fig.[T8]in our effort to reproduce the ID results of[Murphy & Burrows 1 2008) and Nordhaus et aT| <2010D . 

explosion sets in very early (see the black line in Fig. [20b and therefore for large values of the mass accretion rate. The region 
between M « 0.3 M Q s~ 3 and M « 0.8 M s~ 3 is difficult to probe with a stepwise increase of L Ve , because the mass-accretion rate 

there changes so rapidly that the shock shows time-dependent dynamics instead of settling into a quasi-stea dy state. 

Non e of the critical curve s obtain ed with the direct implementation of the neutrino treatment described in IMurphv & Burrowsl 
(120081) and lNordhaus et al.l d2010t) can reproduce the critical luminosity curves reported in these papers reasonably well. We 
therefore decided to proceed with the modifications described in Sect. [2] 
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'pb [s] 

Fig. 20. — Time e volution of the shock radius as a function of post-bounc e tim e, r p i,, for ID simulations performed with the deleptonization scheme of 
Liebendorfei 12005) for electron-fraction trajectory "Deleptonization 1" of Fig. [18] The colors correspond to different electron-neutrino luminosities, which are 
labeled in the plot in units of 10 52 erg s -1 . 
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